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Abstract: 

Coupled triple well {(j)^) one-dimensional potentials occur in both condensed matter physics and field 
theory. Here we provide a set of exact periodic solutions in terms of elliptic functions (domain wall arrays) 
and obtain single domain wall solutions in specific limits. Both topological, nontopological (e.g. some 
pulse-like solutions) as well as mixed domain walls are obtained. We relate these solutions to structural 
phase transitions in materials with polarization, shuffle modes and strain. We calculate the energy and 
the asymptotic interaction between solitons for various solutions. We also consider the discrete analog of 
these coupled models and obtain several single and periodic domain wall exact solutions. 
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1 Introduction 

Many physical problems of interest often have a coupling between two order parameters. In a recent paper 
we obtained exact single and periodic domain wall solutions of coupled double well (0^) models [1, (both 
continuum and discrete) which arise in many second order transitions [21 El El • Similarly, coupled triple 
well or (j)^ models j3] arise in the context of first order phase transitions. There exist analogous coupled 
models in field theoretical contexts (HI |S1. Specifically, when a first order transition is driven by two 
primary order parameters, the free energy should be expanded to sixth order in both order parameters 
with a biquadratic (or possibly other symmetry allowed) coupling. An example of such a transition is 
the cubic to monoclinic structural transformation in shape memory alloys NiTi and AuCd in which two 
phonon (or shuffle) modes are simultaneously active |9|. A simpler example is the square to oblique lattice 
transition driven simultaneously by two different strain components with a biquadratic coupling ^U]. Here 
our motivation is to obtain all possible single and periodic domain wall solutions of these models and then 
connect to experimental observations wherever possible. 

The paper is organized as follows. In Sec. II we provide twenty distinct periodic solutions of a coupled 
continuum (j)^ model. In doing so, we obtain previously unknown periodic solutions of the standard 
(uncoupled) (j)^ model ^2^1^]. We also calculate the total energy and the interaction energy between 
solitons for all the solutions. Section III contains six different solutions for the corresponding discrete cff' 
case. The latter arises in the context of first order phase transitions on a lattice. Finally, we conclude in 
Sec. IV with remarks on related coupled models. 

2 Solutions of a Coupled 0^ Continuum Model 

In the literature, the examples of a coupled (f)^ model are scarce. The reason is presumably that one 
field (or order parameter) is usually dominant and leads to a first order transition. One does not need 
nonlinearity in the secondary order parameter (or the field) to drive the transition. However, in crystalline 
phase transitions this is not necessarily true and within the context of Landau theory, in certain transitions 
symmetry allows one to go to the sixth order nonlinearity in both fields (as primary order parameters) 
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with a biquadratic coupling In this paper, we adopt this viewpoint and consider the solutions of 

a coupled (j)^ model with the potential 

ai ,2 ^1 ,4 , Cl q\ , / «2 , 2 ^2 ,4: C2 A d 2 i2 



v{<t>, i^) = [^-f^' - -f<t>' + -f^"") + [^fr - -fr + -fr) + ^rr , (i) 

where 01.2, ^i,2) ci^2 and d are material (or system) dependent parameters; (j) ^.nd -0 are scalar fields. A 
comprehensive analysis of this coupled continuum (j)^ model shows that it has twenty possible periodic 
solutions, i.e. eleven "bright-bright", three "bright-dark" and six "dark-dark" solutions. In particular, we 
obtain six solutions below the transition temperature Tc, ten at Tc, one above Tc and three in mixed phase 
in the sense that while one of the field is above Tc, the other one is below Tc- The latter situation is akin 
to the one found in multiferroic materials where one transition (e.g. antiferromagnetic) takes place at a 
higher temperature than the other transition (e.g. ferroelectric) or vice versa |141 115j . However, in the 
corresponding discrete case we are able to obtain only six periodic solutions. 
The (static) equations of motion which follow from Eq. are 

d'^d) 
d'^ih 

—r = a2i;- hi^^ + C20^ + dil^cf)'^ . (2) 

These coupled equations have twenty distinct periodic (elliptic function) solutions which we discuss one by 
one systematically. It is worth pointing out that in turn in the single soliton limit, these lead to nine distinct 
(hyperbolic) soliton solutions. In particular, one obtains three solutions below the transition temperature 
Tc, three at Tc, one above Tc and two in the mixed phase in the sense that while one of the field is above 
Tc, the other one is below T^. 

For static solutions the energy is given by 

2 1 / J„;,\ 2 



E 



dx , (3) 



where the limits of integration are from —00 to 00 for hyperbolic solutions (single solitons) on the full line. 
On the other hand, for the periodic solutions (i.e. soliton lattices), the limits are from —K{k) to K{k). 
Here K{k) (and E{k) below) denote complete elliptic integral of the first (and second) kind pT)] . 
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Using the equations of motion, one can show that for all of our solutions 



1 (^\\}^(^\ 



2 \dxj 2 \dx J 



+ C, 



(4) 



where the constant C in general varies from solution to solution. Hence the energy measured with respect 
to the reference energy (e.g. the local or global minimum) E = E — J C dx is given by 



E = E 



Cdx 



dx J \dx 



dx . 



(5) 



Below we give an explicit expression for the energy in the case of all twenty periodic solutions (and hence 
the corresponding nine hyperbolic solutions). In each case we also give an expression for the constant C. 
To begin with, we first obtain ten bright-bright solutions at Tc. 



2.1 Solutions I and II 

We look for the most general solutions to the coupled Eqs. © in terms of Jacobi elliptic functions sn(x, m) 
cn(x,m) and dn(x,m) ^H] where the modulus m = k'^. It is easily shown that 



(j) = A'^ 1 ± sii{Bx + xq, m) , = Fy 1 it sn{Bx + xq, m) , 
are solutions to the coupled Eqs. ((3) provided the following six coupled equations are satisfied 



(6) 



ai - biA^ + dF"^ + ciA^ 

- biA^ + dF^ + 2ciA^ = 
ciA^ 



4 ' 



4 _ 3mB^ 



02 - b2F^ + dA^ + C2F^ 



- b2F^ + dA^ + 2c2F^ 



4 ' 
-niB"^ 



C2F^ 



3mB^ 



(7) 

(8) 
(9) 
(10) 

(11) 
(12) 



Here +(— ) sign in (p goes with +(— ) sign in ip, A and F denote the amplitudes of the kink lattice, B is 
an inverse characteristic length while xq is the (arbitrary) location of the kink. Three of these equations 



determine the three unknowns A, F, B while the other three equations give three constraints between the 
seven parameters ai^2, &i,2) ci^2) In particular, A,F,B are given by 

d2 _ 4ai 2 ^ 8m((i + 62)01 8m{d + bi)ai 

(5m -1)' ~ (5m- l)(6i62 ' ~ (5m - l)(6i62 - d^) ' ^ ^ 

while the three constraints are 

/, , n2 /, , n2 /, , ,2n2 64maici(d + 62)^ /., ,x 

ai = a2, ci d + 62r =C2 d + , {hb2 - = ■ (14) 

3(5m — Ij 

For d = the last constraint reduces to the constraint (for the uncoupled cj)^ [111 1121 113j model) 

bl = . (15) 

^ 3(5m-l) ^ ^ 



So far as we are aware of, the uncoupled (i.e. d = 0) kink lattice solution (p = it sn(i?x + xq, m), 
satisfying 



2 _ 4ai 2 8mai 



(5m - 1) (5m - l)6i ' ^ ' 

and the constraint (|15j) have never been explicitly written down before in the literature. 
At m = 1 these solutions go over to the topological (bright-bright) solutions 



(^ = ^^l±tanh(5x + xo), V = 1 ± tanh(Bx + xq) , (17) 

provided 

2 „ ^2 2{d + b2)ai 2 2{d + bi)ai 
Further, two of the constraints remain the same as above [Eq. (|14|) ] while the third one is now given by 

(6162 - d^f = (16/3)aici((i! + 62)^ , (19) 

which for d = reduces to 

bl = (16/3)aici . (20) 

This corresponds to the point (i.e. temperature) where there are three degenerate minima (or the point 
where a discontinuous transition takes place) of the usual (uncoupled) 0^ field theory jlll 1121 [T^ . 
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Energy: Corresponding to the periodic solutions the energy E and the constant C are given by 



E = ^ — —E{k) , 



C = --{1 - k^){A^ + F^)B'^ . 



(21) 



On using the expansion formulas for E{k) [and K{k) to be used later in the paper] around k = 1 (note 
m = k'^) as given in 



E{k) = 1 + 



k'^ 



In 



k' 



+ 



1 



+ ..., k'^ = l-k'^, 



(22) 
(23) 



.k' J 4 

the energy of the periodic solution can be rewritten as the energy of the corresponding hyperbolic (bright- 
bright) soliton solution [Eq. ((T7)) ] plus the asymptotic (i.e. widely separated solitons) interaction energy 
[ISl. We find 



1 k 



n 



E = Ek^nk + Eint = {A' + F')B - + 



In 



k' 



Note that these solutions exist only when 6162 > d?. The interaction energy vanishes at k"^ 
should. 



(24) 
1, as it 



2.2 Solution III 

It turns out that the coupled Eqs. © also admit a novel mixed solution at T = Tc given by 



AJ 1 + sn(i?x + xq, m) , = FJ 1 — sn{Bx + xq, m) , 



provided the following six coupled equations are satisfied 

ai - biA^ + dF^ + ciA"^ 
- hiA^ - dF^ + 2cxA^ = 



-B' 



02 - b2F' + dA^ + C2F^ 



b2F^ - dA^ + 2c2F^ 



4 ' 
-mB'^ 



-B' 

4 ' 
-niB'^ 



(25) 



(26) 

(27) 
(28) 
(29) 
(30) 



c.F* = ^. (31) 

Three of these equations determine the three unknowns A,F,B while the other three equations give three 
constraints between the seven parameters 0-1^2^^1,2, ci^2jd. In particular, A,F,B are given by 

^2 4ai(6i&2-a!^) 



[(5m - 1)6162 - IGmbid + (11m + l)^^)] ' 

^2 ^ 8mai(62 - d) ^ 8mai(6i - d) 

[(5m - 1)6162 - 16mbid + (11m + l)d'^)] ' [(5m - 1)6162 - 16m6id + (11m + 1)^2)^ ^ 

while the three constraints are 

[(5m - 1)6162 + (11m + l)d'^ - 16mb2d]ai = [(5m - 1)6162 + (11m + l)d'^ - 16mbid]a2 , 

ih (u tu u 64maici(62 - d)^ , . 

ci(62 -d) = C2(6i - d) , (6162 -d) = 3[(5^_i),^,^^(nm + l)d2-16m6id] " ^'^^ 

For d = the last constraint reduces to the constraint (for the uncoupled (j)^ model) as given by Eq. (|15() . 
Note that no new solution is obtained by having the — sn term inside the square root in cj) and the +sn 
term inside the square root in ip in Eq. H25|) since such a solution is trivially obtained from Eq. (|25() by 
merely interchanging A and F. This is so because the potential V{(j), ip) is completely symmetric in (j) and 
ip. Similar comments apply to many of the other asymmetric solutions given below. 
At m = 1 this solution goes over to the topological (bright-bright) solution 



4) = A^l + tanh(5x + xq) , il) = F^l- tanh(5x + xq) , (34) 

provided 

2 ^ 01(6162 - d^) 2 ^ ai{b2-d) 2 ^ - d) 

(6i62-6i(i + d2) ' 2(6162 -46id + 3^2) ' 2(6162 - 46id + 3^2) ' ^ ' 

while the three constraints are 

(6162 + - Ab2d)ai = (6162 + 2>d^ - Abid)a2 , 

can - ^r- ^ - . (M. - <e) ^ 3(m:':'3^-1Co ■ '^^> 

For d = 0, the last constraint reduces to the constraint ()2U|) which is the point where there are three 
degenerate minima (i.e. the point where a discontinuous transition takes place) of the usual (uncoupled) 
06 field theory [TTIIIIIISI- 



Energy: Remarkably, corresponding to the periodic solution ()25() . the energy E is again the same as 
that for solutions I and II and is given by Eq. Thus the interaction energy too remains unchanged 

and is again given by Eq. (|24|) . However, C is different for this solution than those for solutions I and II 
and is given by 

C=<^i^^±^-|(^= + F=)(3 + ™). (37) 

Note that this solution also exists only when 6162 > d?- 

2.3 Solutions IV and V 

It is easily shown that 



(j) = Ay 1 lb ^/rnsn{Bx + xq, m) , = Fy 1 it ^/rnsn{Bx + xq, m) , (38) 
are also solutions to the coupled Eqs. ^ provided the following six coupled equations are satisfied 

ai - biA^ + dF^ + ciA^ = , (39) 

- hiA^ + dF^ + 2ciA^ = ^ , (40) 

ciA' = ^ , (41) 

a, - b^F^ + dA^ + C2F^ = , (42) 

- b2F^ + dA^ + 2c2F^ = , (43) 

C2F' = ^ . (44) 

Three of these equations determine the three unknowns A, F, B while the other three equations give three 
constraints between the seven parameters ai^2-,bi,2-,ci^2-,d. In particular, A,F,B are given by 

2 ^ 4ai 2 ^ 8(^ + 62)01 p2 ^ 8{d + bi)ai 

(5-m)' (5-m)(6i62 -d2) ' (5 - m)(6i62 - d^) ' ^ 

while the three constraints are 

a, = a2, c^{d + b2? = C2{d + h)\ ihh2-d'f = ^-^^^^±^. (46) 



For d = the last one reduces to the constraint (for the uncoupled (p^ model) 1121 113j 



So far as we are aware of, the uncoupled (i.e. d = 0) kink lattice solution (f) = Ay 1 it ^/rnsn{Bx + xq, m), 
satisfying 

(5 - m) ' (5 - m)bi ' ^ ^ 

and the constraint ()47() have never been explicitly written down before in the literature. 

At m = 1 these solutions go over to the topological (bright-bright) solutions H17|) satisfying the con- 
straints (HTni to CT. 

Energy: Corresponding to the periodic solutions H38|) . the energy E and the constant C are given by 

C = ^{l-e){A'^ + F^)B\ (49) 

On using the expansion formulas for E{k) and K(k) around k = 1 (note in = k'^) as given by Eqs. (|22() and 
H23|) . the energy of the periodic solutions can be rewritten as the energy of the corresponding hyperbolic 
(bright-bright) soliton solutions [Eq. lfT7|) ] plus the interaction energy. We find 



E = Ekink + Eint = {A + F )i? ( — 



(50) 

Note that these solutions also exist only when 6162 > d^ . The interaction energy vanishes at /c = 1. 



2.4 Solution VI 

It turns out that the coupled Eqs. ((2j also admit a novel mixed solution at T = Tc given by 



cj) = 1 + ^/msxl{Bx + xo, m) , = -^y 1 ~ ^/ms\i{Bx + xq, m) , (51) 
provided the following six coupled equations are satisfied 

ai - 61^2 + ^p2 ^ ^^^4 ^ , (52) 
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ciA' = ^ , (54) 

a2 - 62^2 + dA^ + C2F^ = , (55) 

- 62^2 - dA^ + 2c2F^ = ^ , (56) 

= ^ . (57) 

Three of these equations determine the three unknowns A,F,B while the other three equations give three 
constraints between the seven parameters ai^2,bi.2,ci^2,d. In particular, A,F,B are given by 



,2 4ai(6i62-d^) 



[(5 - m)bib2 - Wbid + (11 + m)d^)] ' 

.2 ^ 8mai{b2 - d) 2 ^ 8mai(6i - d) 

[(5-m)6i62-166id+(ll + m)d2)] ' [(5 - m)6i62 - 166id + (11 + m)d2)] ' ^ ' 

while the three constraints are 

[(5 - m)bih2 + (11 + m)d^ - lQb2d\ai = [(5 - m)bib2 + (11 + m)d^ - 166id]a2 , 

For (i = the last one reduces to the constraint 1)47^ . 

At m = 1 this solution goes over to the topological (bright-bright) solution H34|) satisfying the same 
constraints as given by Eqs. and 

Energy: Remarkably, corresponding to the periodic solution (|51() . the energy is again the same as 
that for solutions IV and V and is given by Eq. H49() . Thus the interaction energy too remains unchanged 
and is again given by Eq. 1(50]) . However, C is different for this solution than those for solutions IV and V 
and is given by 

c = — - T6 + ^^^^ + ^""^ ■ 

Note that this solution also exists only when 6162 > cP ■ 
2.5 Solutions VII and VIII 

Apart from the above six solutions, at T = Tc, the field equations also admit mixed type of solutions. In 
particular, it is easily shown that 



A\J \ rtsn{Bx + XQ^rn) , il) = Fy 1 ± y/rnsn{Bx + XQ,m) , (61) 
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are solutions to the coupled Eqs. ^ provided the following six coupled equations are satisfied 

_ 732 

ai - biA'^ + dF^ + ciA^ = , (62) 
- hiA^ + ^dF^ + 2ci^4 = , (63) 

ci^' = ^^, (64) 
as - 62^2 + + ^2^4 ^ , (65) 

- ^^2^' + ^i^' + 2V^C2F4 = ^^"^ , (66) 
C2F^ = ^ . (67) 

Three of these equations determine the three unknowns A, F, B while the other three equations give three 
constraints between the seven parameters ai^2-,i'i,2^ci,2^d. In particular, A,F,B are given by 

^2^ 4ai(6i62 -d^) 

A^ = 
F^ = 

while the three constraints are 

[(5m - 1)6162 - (8^/^ - 3m - l)d'^ - 8dbi{l - Vm)]a2 

= [(5 — m)6i62 — (8\/m — 3 — m)d'^ + 8db2V^{l — \/rn)]ai , ci(d + \/m^2)^ = £2(61 + \/rrid)'^ , 

(h h - d^) = 64aici(d + Vm62)^ , , 

' 3[(5m - 1)6162 -(8V^- 3m - l)d2 - 8d6i(l - V^)] ■ ^ ^ 

For d = the last constraint reduces to the constraint ()15() . 

At m = 1 these solutions go over to the topological (bright-bright) solutions ((T7|) satisfying the con- 
straints (HTni to CT. 

Energy: Corresponding to the periodic solutions (|61|) . the energy £^ and the constant C are given by 

C = 1 [ai^2 ^ ^^^2j _ :^ [(3 + A:2)^2 ^ ^ 3^2)^2^ _ (70) 
4 16 
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[(5m — 


1)6162 - ( 

< 


8y/m — 3m — l)d'^ — 
i^/m{d + ^/mb2)ai 


8d6i(l - ■ 






[(5m — 


1)6162 - ( 


8^/m — 3m — 1)^2 — 


8dbi{l - , 










8(61 -|- y/md)ai 






(68) 


[(5m — 


1)6162 - ( 


8y/m — 3m — l)d'^ — 


8d6i(l - ■ 





On using the expansion formulas (|22() and ()23() . the energy of the periodic solutions can be rewritten as 
the energy of the corresponding hyperbolic (bright-bright) soliton solutions [Eq. (|17jl ] plus the interaction 
energy. We find 



E = Ekink + Eint = {A^ + F^)B ( ^ + V 



1 k'^ [ /4\ 1 

^"'Fj-2 



BF'^ , ,o / 4 



-k"ln[-]. (71) 



Note that these solutions exist only when 6162 > d?- The interaction energy vanishes at fc = 1, as it should. 
2.6 Solutions IX and X 

Finally, there are two more mixed type of solutions at T = given by 



(j) = Ay 1 lb sn{Bx + xq, m) , t/j = Fy 1 =F y/rnsn{Bx + xq, m) , (72) 
provided the following six coupled equations are satisfied 

_ 732 

ai - biA^ + + ciA^ = , (73) 
- hiA^ - ^dF^ + 2ci^4 = , (74) 

cU^^^, (75) 
a2 - 62^2 + dA^ + C2F4 = , (76) 

- ^^62^2 - + 2V^C2F^ = , (77) 

C2F^ = ^ . (78) 

Three of these equations determine the three unknowns A, F, B while the other three equations give three 
constraints between the seven parameters ai,2,bi,2,ci^2,d. In particular, A,F,B are given by 

^2 _ Aai{bib2 - d"^) 

A^ 



[(5m 


- 1)6162 + (11m + l)c?2 - 16^ 
8^/m{^/mb2 — d)ai 


frndbi] 


[(5m 


- 1)6162 + (11m + 1)^2 -16^ 
8y^(6i — ^/md)al 


/mdbi] 



[(5m - 1)6162 + (11m + 1)^2 - 16Vmd6i] ' ^ ^ 
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while the three constraints are 

[(5m - 1)6162 + (11m + - 16d^yTEbi]a2 = 
[(8-v/m — 3 — 771)6162 + (8-v/m + 3 + m)cf — 16md62]oi , 
Ci{\/rnb2 — d)"^ = mc2(6i — ^/mdf , (6162 — d^) 



2_ ^j\2 /. > ,2n _ 64aici(^62 - d)^ 



3[(5m - 1)6162 + (11m + l)d^ - Wd^h] ' 

(80) 

For d = the last constraint reduces to the constraint (|15() . 

At m = 1 these solutions go over to the topological (bright-bright) solutions (pi)) satisfying the con- 
straints and ((Snj- 

Energy: The energy E and the constant C for the solutions (|72() are the same as for solutions VII and 
VIII and are given by Eq. H70() and hence the interaction energy is again as given by Eq. ()71|). Note that 
these solutions also exist only when 6162 > d^. 

2.7 Solution XI 

We shall now present six periodic solutions which exist in the case T < Tc- Out of these one is bright-bright, 
three dark-dark and two bright-dark solutions. It is easily shown that 

^ Asn{Bx + xo,m) Fsn{Bx + xo,m) 



\/l — Dsn^{Bx + XQ,'m) ' ^yl — Dsn^{Bx + xq, m) 

is a periodic variant of the bright-bright solution provided the following six equations are satisfied 

ai = (3D - 1 - m)B'^ , (82) 

dF'^ - biA^ - 2aiD = 2(m -D- mD)B'^ , (83) 

- dDF"^ + hiDA^ + aiD^ + ciA^ = mDB'^ , (84) 

02 = (3D - 1 - m)B^ , (85) 

dA^ - 62^^ - 2a2-D = 2(m - D- mD)B'^ , (86) 

- dDA^ + 62DF2 + aaD^ + 02^^ = mDB^ . (87) 
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Prom here one can determine the four unknowns A, F, B, D and further one obtains two constraints between 
the seven parameters ai_2-&i,2) ci_2, d. We obtain 

„2 ai ^2 2{d + h2)B'^[2{l + m)D -m-2,D'^] 

= , r , A — 



F 



(SL* - 1 - m) ' (6162 - <P) 

2 _2{d + bi)B'^[2{l + m)D - m - 3D^] 
~ (6162 - d^) ' 

[2(1 +m)D - m-3D^]^ 3(6162 - d^)^ 



D{1- D){m- D){3D -1-m) 4aiCi(d + 62)^ ' 
while the two constraints are 

ai=a2, C2{d + hif = ci{d + h2f . (89) 

Using > as well as D < m, it follows that 

1 + m — \/l — m + m? 



<D <m. (90) 



As a result, ai, 02 are positive so long as 



1 + m — — m + m? ^ 1 + m , . 
7, < D < — - — , (91) 



while they are negative if 

1 + m 



<D<m. (92) 



3 

In the limit d = 0, one obtains the uncoupled kink lattice solution jlll 1121 [T^ 

^ Asn{Bx + XQ,m) 

Y/r^-"Dsn^(^x~+~xo7"i) 

satisfying 

2 2B^[2D{l + m)-m-3D^] 2 «i 



A^ = ' , , B' 



hi ' (3L> - 1 - m) ' 

[2D(1 + m)-m- 3D'^]'^ _ 35f 



D{1 - D){m - D){3D - 1 - m) 4aiCi 
At m = 1 the solution reduces to the topological (bright-bright) solution 



1 - i:>tanh''(Bx + xo) a/1 - tanh''(-Bx + xq 



(93) 



(94) 



Atanh(Bx + xq) F tanh.(Bx + xq) 

, V = ; , (95) 
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provided 



2 ^ Qi ,2 ^ 2{d + b2)B^{l-D){3D-l) 

{3D -2)' {bih2-d?) 
^2 2{d + hi)B'^{l- D){3D -I) 



(6162 - 



3 



1 + 



(6162 - f^' 



'2\ 



(96) 



V(&ife2-rf')'-4aiCi(d + 62)^- 
while the two constraints are the same as given by Eq. H89|) . Since > 0, it then follows that 1/3 < D <1. 

Further, ai,a2>0if2/3<Z)<l while 01,02 < in case 1/3 < D < 2/3. The constraint D <1 implies 

the inequality 

(6162 - cff > (16/3)aici((i + 62)^ , (97) 

which for d = reduces to 

bl > (16/3)aici . (98) 

This corresponds to the situation when there are two degenerate absolute minima and a local minimum or 
maximum at </> = depending on if oi > or oi < 0, respectively [TTl [T^ [T^ . 
Special case of 6162 = d'^ 

One can show that the solution (|81() exists even in the case 6162 = d^. It turns out that such a solution 
exists only if 

biA^ = b2F^ , (99) 

and further 

,2 _ -oi 



oi = 02 < , 61 C2 = 62C1 , B 



\/l — m + m 



2 ' 



D = ^— — ^ ^- , ciA"^ = B^{1- D)[m{l + D) -2D]. (100) 

In the limit m — > 1, relation (|100j) takes the simpler form 

B' = -au D = l, c^A' = -^. (101) 

Energy: Corresponding to the periodic solution H81|) . the energy E and the constant C (using the 
appropriate integrals given in ^7j) given by 

E = {A'^ + F'^)BIi , 
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where 



with 



C = -t{A^ + F^)B^ 



[AiKjk) + A2E{k) + A^niP, k)] 
4D^{1- D){k^ - D) 



(102) 



(103) 



Ai = {k^ - D){2D + k^ - 3D^) , A2 = D[2D{1 + k^) - k^ - 3D^] , 

As = 4(1 - D){k'^ - D)[{1 + k^)D - k^] - [D^ - 2(1 + k'^)D + 3k^][2{l + k^)D - k^ - 3D^] (104) 

Here 11(1), k) denotes the complete ehiptic integral of the third kind |161 117j . 

In order to obtain the interaction energy corresponding to this solution, one needs the expansion formula 
for n(D, k) around k = 1. Unfortunately, only the leading term in the expansion is given in various books 
that we have encountered, but the subleading term is not mentioned. In particular, the result well known 
in the literature is 



a 



K{k = 1) - (1 - a^)n(a^ A; = 1) = - In 



1 + a 
1 — a 



(105) 



We now show that starting from the basic definitions of K{k) and n(a^,/c), the subleading term in the 
expansion is easily derived. 

Around A; = 1, we have (note k'"^ = \ — k"^) 



K(k) - (1 - a2)^(a^ k) = K{1) - (1 - a^)n{a^, 1) + {k^ - 1) 
where K{k) and n(a^,A:) are defined by ITT) 



K{k) 



dx 



dk^ ^ " 'dk^ 



dx 



fc2 = l 



+ 0(A:'^), (106) 



(107) 



/o y/{l-x')il-k''x') ' ' ' ' 7o {l-a^x^)^{l-x'){l-k'x') ■ 
The well known result (|1U5|1 is easily obtained by using basic definitions (|107jl of K(k) and 11(0^, A;), i.e 



K(l)-(l-a2)^(a^l) =q2 _ 



dx 



a 



In 



1 + a 



1 — a 



(108) 



Similarly, using the basic definitions (|107)) it is straightforward to show that 



dK 
dk^ 



dli 

'dA^ 



a 



dx 



dx 



fc2=i 2(1 - a2) Uo (1 - x2) 7o (1 - a2x2^ 
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(109) 



These integrals are easily evaluated and thus we find that 



a 



K{k) - (1 - a^)U{a, k) = - In 



1 + a 



l-a 



2(1 -q2) 



In 


4' 


-im 


'1 + a' 








) 




k'. 


2 


1 — a_ 





+ 0{k'^) 



(110) 



On using the expansions of E{k), K{k), and Il(D,k) around k = 1 as given by Eqs. H22|) . I\2'6\i and 
H11U|) . the energy of the periodic solution can be rewritten as the energy of the corresponding hyperbolic 
(bright-bright) soliton solution [Eq. (|95j) ] plus the interaction energy. We find 



E = Ek,nk + E,nt = {A^ + F^)B[lf'^ + k'^l[^^] , 



(111) 



where ij*^^ and are given by 



r(0) 



{3D-1) ^(l^aD)j^ 



AD{l-D) 81)3/2 



I 



(1) 



:i+3L>2) 



1 



In 



1 + Vd' 

1-VD 

'i + Vd 



1-^ 



(112) 
(113) 



1 16D{1 - Df 161)3/2 
Note that this solution exists only when 6162 > ■ The interaction energy, as expected, vanishes at 
k = l. 



2.8 Solution XII 

There are actually three periodic solutions, all of which at m = 1 reduce to the same hyperbolic nontopo- 
logical (dark-dark) soliton solution. However, for m < 1, all three (periodic) solutions are distinct which 
we now discuss one by one. One of the solution is 

Ac\'i{Bx + XQ^m) ^ Fcn{Bx + XQ,m) 



\/l — Dsn'^{Bx + xo,m] 
provided the following six equations are satisfied 



^1 -Dsn2(Bx + 



xo,m) 



ai - biA^ + dF^ + ciA'^ = -(1 - D)B'^ , 



(1 + D)biA^ - (1 + D)dF^ - 2aiD - 2ciA^ = 2(m - D){1 - D)B^ 
dDF'^ - biDA^ + aiD"^ + ciA^ = mD{l - D)B^ , 
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(114) 



(115) 

(116) 
(117) 



a2-b2F^ + dA^ + C2F^ = -{I- 0)3"^ , (118) 

(1 + D)b2F'^ - (1 + D)dA'^ - 2a2D - 2c2F^ = 2{m - D){1- 0)3"^ , (119) 

dDA^ - b2DF^ + a2D^ + C2F^ = mD{l - D)B'^ . (120) 

From here one can determine the four unknowns A, F, B, D and further one obtains two constraints between 
the seven parameters 01^2, ^1,2, ci^2, d. We get 

.2 2{d + b2)B^[{l + 2D)m- {D + 2)D] ^2 2{d + bi)B^[2{l + 2D)m - {D + 2)D] 

^ = n T^\n , W\ ' ^ ~ 



[l-D){b,b2-d^) ' (l-D)ibib2-d^) 

^2 ^ il-D)ai [jl + 2D)m- (2 + D)D]^ ^ 3(6ib2 - d^f 

[2m - 1 - Z)(2 - m)] ' D{m - D)[2m - I - D{2 - m)] 4aiCi(d + 62)^ ' 

while the two constraints are 

ai=a2, C2[d + bif = ci{d + b2f . (122) 

Using > 0, it fohows that 



Q < D < y/l - m + m? - {I - m) . (123) 

It is worth noting that for m < 1/2, oi, 02 < irrespective of the value of D. 
In the hmit d = 0, one obtains the uncoupled pulse lattice solution 

Acn{Bx + XQjTn) 



satisfying 



\/l — Dsn'^{Bx + xo,m) 



2B^[{l + 2D)m- {D + 2)D] o (1-15)01 



(124) 



{1-D)bi [2m - 1 - D{2 - m)] 

[{l + 2D)m- {2 + D)D]'^ 35f 



(125) 



D{m - D)[2m - 1 - D{2 - m)] 4aici 

So far as we are aware of, this solution has never been explicitly written down before in the literature. 

In the limit m = 1, the solution H114|l goes over to the hyperbolic, nontopological, dark-dark soliton 
solution 

^sech(i3x + xo) , Fsech(i?a; + xo) 

:, V' = ; (126) 



1- Dtanh-^iBx + xo) Jl - Dtanh^{Bx + xq 
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provided 

{hih-d?) ' D 4aiCi(d + 62)2 ' ^ ' 

while the two constraints are the same as given by Eq. (|122() . Note that ai,a2 are now positive. The 
condition D <1 gives the same constraint as given by Eq. ()97|) . 
Special case of 6162 = (f 

One can show that the solution (|114|) exists even in the case 6162 = d^- It turns out that such a solution 
exists only if 

biA" = b2F^ , (128) 

and further if 



ai = 02 < , 61C2 = 62C1 , B'^ — 



\/l — m + m 



2 ' 



D = \/l — m + rn^ — (1 — m) , ci^^ = — ai 



(2 - m) 



y/l — m + m? 



(129) 



In the limit m — > 1, the constraint (|129() gives 

B"^ = -ai, D = l, ci^^ = -ai . (130) 
Thus, for the case 6162 = the hyperbolic dark-dark soliton solution (|126|) reduces to a constant solution 

(l) = A, i> = F. (131) 

Energy: Corresponding to the periodic solution (|114j) . the energy E and the constant C are given by 

E = {A^ + F'^)Bl2 , 

2(1 -D) ^ ^ 



Here 



, (1 - D)[BiKik) + B^Ejk) + B^UjP, k)] 
AD^k^-D) ' 
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with 



-(k^ - D)[(l - D)^ - [1 - k^){l - AD)] D[{l-D'-{l-k^){l + 2D)] 

0^:0^ ' 0^^^ 

= 4(.^ - Dm + D)k^ - D] - \£l^li^^ + - - - D^] 



(134) 



(l-D) 

On using the expansion formulas for E{k), K{k) and 11(1), k) around = 1 as given above, the energy of 
the periodic solution can be rewritten as the energy of the corresponding hyperbolic (dark-dark) soliton 
solution [Eq. (|126|) ] plus the interaction energy. We find 



E = Ek,nk + Eint = {A^ + F^)B[l'i^ + k'^li^^] , 



(135) 



where l2^^ and Iq^^ are given by 



,,0,_(1 + O)_(lz£)!,„ 



4D 



81)3/2 



1 + Vd' 
1 - Vd 



(1) ^ (Z)2 -AD-1) (31)2 + QD_i) 
2 16D{l-D) 16Dy^l-D) 



1 + Vd 
1 - Vd 



1 , /4 
+ ^In 



[l-D] \k' 



(136) 



(137) 



Note that this solution exists only when 6162 ^ d?- The interaction energy, as expected, vanishes at 
k = l. 



2.9 Solution XIII 



Another solution is given by 



Adn(Bx + XQ,m) 

(f) = — 

yjl — DsD?{Bx + Xq, m) 

provided the following six equations are satisfied 



Fdn[Bx + xo, m) 
yj\ — Ds\i^{Bx + xo, m) 



(138) 



ai - hiA^ + dF'^ + ci^^ = -(m - D)B'^ , 

(m + D)biA^ - (m + Z))^^^ - 2aiD - 2mciA'^ = 2{m - D){1 - D)B^ , 
dniDF"^ - biniDA^ + aiD^ + cim^yl^ = D{m - D)B'^ , 
a2 - b2F^ + _^ g^F^ = -(m - i:))^^ , 



(139) 

(140) 
(141) 
(142) 
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(m + D)b2F'^ - (m + D)dA'^ - 2a2D - 2mc2F^ = 2{m - D){1 - D)B'^ , (143) 
dmDA^ - b2mDF'^ + asD^ + cgm^F^ = D{m - D)B'^ . (144) 

From here one can determine the four unknowns A,F,B,D and further one obtains two constraints between 

the seven parameters ai_2, &i,2) ci^2, d. We obtain 

2 _ 2(d + b2)B^[2D{l - m) + m - L>2] ^ ^ 2{d + bi)B^[2D{l - m) + m - D^] 

~ {m - D){bib2 - d'^) ' " {m - D){bib2 - d^) ' 

2 ^ (m - £')ai [2D{l-m) + m-D^]^ ^ 8(6162 - d^)^ 

[2m(l - D) - m2 + ' L>(l-i:))[2m(l-D)-m2 + Z)] 4aiCi((i + 62)2 ' ^ ^ 

while the two constraints are 

ai = a2, C2[d + bif = ci{d + b2f . (146) 

The best bound on D is given by < D < m which ensures > 0. Further, unUke two of the above 
solutions (i.e. solutions XI and XII), ai,a2 are always positive. 

In the limit d = 0, one obtains the uncoupled pulse lattice solution 

Adn{Bx + xo,m) 



satisfying 



^1 — Dsn^{Bx + xq, m) ' 

2B^[2D{l-m) + m- D"^] {m - D)ai 

B — 



(147) 



{m-D)bi ' [2m{l - D) + D - m?] 

[2D{l-m)+m- D'^f 36? 



(148) 



D{1 - D) [2m{l -D) + D-m?] AaiCi ' 
In the limit m = 1, the solution H138|) also goes over to the hyperbolic, nontopological, dark-dark soliton 
solution (|126|) satisfying the same constraint as given by Eq. (|127jl . 
Special case of 6162 = 

One can show that the solution (|138|) exists even in the case 6162 = d^. It turns out that such a solution 
exists only if 



and further if 



2 

ai = a2 , 61C2 = 62C1 , B — 



\/l — m + m? 



D = ^l-m + m? + il-m), ciA^ = — . (150) 

(m — D) 
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In the limit m ^ 1, the constraint (|15(J() . as expected, reduces to Eq. (|130() and the solution (|138|) reduces 
to the constant solution (|131|) . 

Energy: Corresponding to the periodic solution (|138j) . the energy E and the constant C are given by 

2(A;2 - D) • ^^^^^ 



Here 



with 



(tf - D)[DiK{k) + P,E(t) + P3n(P.fc)] 



p, = _ uf - (1 - . p, = Mlz + (1 - '^)(^" - 1)1 



Pa = 4(1 - P)|(l - P).-' + P| - - ^" + ^''^ + ?f + ^'^ - ^'>° - . (153) 



(A;2 _ 

On using the expansion formulas for E{k), K{k) and 11(1?, A;) around k = 1 obtained above, the energy of 
the periodic solution can be rewritten as the energy of the corresponding hyperbolic (dark-dark) soliton 
solution [Eq. (|126|) ] plus the interaction energy. We find 

E = Ek^nk + Eint = {A^ + F^)B[4^'> + k'^4^^] , (154) 
where I^^ is in fact the same as l2'\ and given by Eq. (|136|) while /g^^ is given by 



(1) _ (L)2 + _ 5) (71^2 _ 2D + 3) 
3 16D{l-D) 16L>3/2(i_£)) 



1 + Vd 
1 - Vd 



Infi). (155) 



[1-D) \k 



Note that this solution exists only when 6162 > cP- As expected, the interaction energy vanishes at 
k = l. 

2.10 Solution XIV 

Yet another solution is 

Acn.{Bx + xo,m) F(hi{Bx + XQ,m) 



a/1 — Dsn^iBx + xq, m) ' ^1 — Dsn-^{Bx + xq, ?n) 
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4ci [{m - 


-D^){d + h2)-2D{l- 


m){b2 - 




2^2 [(m 


-D^){d + b2)-2D{l- 


- m){b2 


-d)] 




{bih2-d?)[l-D) 






2B\m 


- D^){d + hl) + 2D{l- 


- m){bi 


-d)] 



provided the following six equations are satisfied 

ai-biA^ + dF'^+ciA'^ = -{l-D)B^, (157) 

(1 + D)biA'^ - (m + D)dF'^ - 2aiD - 2cxA^ = 2(m - D){\ - D)B^ , (158) 

dmDF^ - biDA^ + aiL>2 + ciA^ = mD{l - D)B'^ , (159) 

a2-b2F'^ + dA'^ + C2F^ = -{m- D)B^ , (160) 

(m + D)b2F^ - (1 + D)dA^ - 2a2D - 2mc2F^ = 2{m - D){1 - D)B'^ , (161) 

dDA^ - b2mDF^ + aaZ?^ + C2m^F^ = B{m - D)B'^ . (162) 

Solving these equations, we determine the four unknowns A, F, B, D 

^2 3D{1 - D){m - D){bib2 - d^f 

A' 

^ ^ {bib2 - d'^){m - D) 

(1 - D)ai + (1 - m)dF'^ = [2m - 1 - (2 - m)D\B'^ , (163) 

while the two constraints between the seven parameters are 

(1 - DfciA"^ = {m- Dfc2F^ , (164) 

and 

1 - D)a2 = 1 - m)b2F^ + ^ ^— ^— . 165 

(m — U) 

In this case too the best bound on D is given by < < m. Further, ai < in case m < 1/2. 

In the limit m = 1, the solution H156|) also goes over to the hyperbolic, nontopological, dark-dark soliton 
solution H126|) satisfying the same constraints as given by Eq. (|127() . 

Unlike the solutions XI, XII and XIII, the solution H156|) does not exist in the case 6162 = d?- This is 
because, in order that such a solution exist, m and D must satisfy 

V^[m - _ 2(1 _ m)D] = -v^[m - D"^ + 2(1 - m)D] , (166) 
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which is impossible. 

Energy: Corresponding to the periodic solution ()156jl . the energy E and the constant C are given by 



E = B{A'h + F'h) , 

b 



c ' 



2 3 2 3 



167) 



2^3 

where I2,l3 are as given by Eqs. ()133() . H134|) . (|152|) and 1)153^ . On using the expansion formulas for 
E(k), K{k) and n(D, k) around A; = 1 as derived above, the energy of the periodic solution can be rewritten 
as the energy of the corresponding hyperbolic (dark-dark) soliton solution [Eq. (|126|) ] plus the interaction 
energy. We find 

E = Ekink + Eint = [{A^ + F^)!^^ + + F^lll^B , (168) 

where I2^\l2^\l^^^ are given by Eqs. ()136() . H137() and l|155() respectively. 

Note that this solution exists only when 6162 > cP- As expected, the interaction energy vanishes at 
k = l. 

2.11 Solution XV 

Finally, we obtain two periodic solutions which at m = 1 go over to the same bright-dark soliton solution. 
The first solution is given by 

^_ Asn{Bx + xo,m) ^_ Fcn{Bx + xo,m) {169) 

y^l — Dsii^{Bx + xo, m) y^l — Ds\i-^{Bx + xq, m) 

provided the following six equations are satisfied 

ai + dF'^ = {'iD-l- m)B'^ , (170) 

biA^ + d{l + D)F'^ + 2aiD = 2{D + Dm - m)B'^ , (171) 

dDF^ + biDA^ + aiD^ + ci^^ = mDB^ , (172) 

02 - 62F2 + C2F^ = -{l-D)B^, (173) 

(1 + D)b2F'^ + dA^ - 2a2D - 2c2F^ = 2(m -£>)(!- D)B^ , (174) 



24 



- dDA'^ - b2DF'^ + a2D'^ + C2F^ = mD{l - 0)3"^ 



(175) 



On solving these equations we determine the four parameters A,F,B,D as well as obtain two constraints 
between the seven parameters. We find 



A 



2 _ 2B'^[b2{2mD -m + 2D- 3D^) - d{2mD + m - 2D - D^)] 

{bib2 - (P) 

^2 _ 2B'^[bi{2mD + m - 2D - D^) - d{2mD - m + 2D - SZ)^)] 
" {1 - D){blb2 - d^) 

2_ 3D{1 - D){m - D){bib2 - d^)^ 



(1 + m) — \/l — m + m? i — 
< D < \J\ 



'1 -Z)tanh^(Bx + 2;o) y^l - Z)tanh^(5x + xq 

and the relations ()176|) to ()179() take the simpler form 

,2 2i?2(l-i?)[62(3Z)-l)-d(l + Z))] 



A' 



(6i&2 - d^) 



(176) 



~ 4:Ci[b2{2mD - m + 2D - SD"^) - d{2mD + m - 2D - D^)]'^ ' 
ai + dF^ = [3i:> - 1 - m]B^ , (177) 

while the two constraints are 

ciA* = (I - Dfc2F\ (178) 
„,-,,f^^ «^l'^-°);_+^^('"-°)l ,0. (179) 

From positivity considerations, one can show from here that 



_ . _ m + — (1 — m) . (180) 

3 

Further, it follows that oi < if I) < (1 + m)/3. 

At m=l, the solution H169|) goes over to the hyperbolic, bright-dark soliton solution 

^ Atanh{Bx + xq) Fsech{Bx + xo) nsT\ 

1 V / • l-'-^-'-j 



j,2 2BHl-D)[bi{l + D)-d{3D-l)] 



^2 _ mbib2 - dY 



AciMSD - 1) - d{l + D)]^ ' 
ai + dF^ = {3D - 2)^2 ^ (ig3) 
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while the two constraints are 

ciA^ = (1 - Dfc2F^ , (184) 

02-62^2 + 52(1 + 21?) = 0. (185) 

Prom here one can show that 1/3 < -D < 1 and further ai <0ifI?>2/3 while 02 < if D > ^ • 
Special case of 6162 = d? 

One can show that the solution H169|l exists even in the case 6162 = d^. It turns out that such a solution 
exists only if 

a/6^[2(1 + m)D -m- SD^] = ^^[m - - 2(1 - ■m)D\ , (186) 

and further while relations (|177|) to (|179|) are still valid, the relations l|176|) are no longer valid. Instead of 
separate expressions for A^, F"^, one now only has a constraint 



biA^ + VM2(1 - D)F^ = 2B^[2{1 + m)D - 3D^ - m] . (187) 

In the limit m — > 1, the hyperbolic bright-dark solution 1)1811) is still valid provided D has the value 

D = , (188) 

3v62 - V Oi 



while the constraint p87() takes the simpler form 

biA^ + \/M2(1 - D)F'^ = 25^(1 - D){'iD - 1) . (189) 

The relations (|TH5|) to (fTH5|) are still vahd. 

Energy: Corresponding to the periodic solution ()169() . the energy E and the constant C are given by 



E = B{A^h + F^h) , 

h-. ^ . _ „ _ A„ _ ^„ . 1 

(190) 



c ' 



A\aiD' + ^DA^ + ^A^ + dDF^) + F\-a2D^ + ^DF^ - ^F^] 
^ 2 3 ^ ^ ^ 2 3 ^ 



21)3 

where Ii,/2 are as given by Eqs. (|1U3() . (|1U4() . ()133|) and (|134|) . On using the expansion formulas for 
E{k), K{k) and n(D, k) around = 1 as derived above, the energy of the periodic solution can be rewritten 
as the energy of the corresponding hyperbolic (bright-dark) soliton solution [Eq. (|181|) ] plus the interaction 
energy. We find 

E = Eunk + E^nt = BA^[l[''^ + A:'2/f )] + iJF^f/f + fe'^jW] ^ (191) 
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where /{J^^ are given by Eqs. TTHh . (tM . (ITO and (tHTI) . 

Note that this solution exists only when 6162 > (i^- As usual, the interaction energy vanishes at /c = 1. 

2.12 Solution XVI 

Finally, the last solution below Tc is given by 

^_ Asn{Bx + xo,m) Fdn{Bx + XQ,m) , , 



^/T^^lysI?{W!F+lcQ~m) ' \/l — Dsn'^{Bx + xq, m) ' 

provided the following six equations are satisfied 

ai + dF^ = (3i:> - 1 - m)B'^ , (193) 

61^2 + d{m + D)F'^ + 2aiL> = 2{D + L>m - m)B'^ , (194) 

dmDF^ + biDA'^ + aiD^ + ci^^ = mDB'^ , (195) 

a2-b2F^ + = -{m- D)B'^ , (196) 

(m + D)62i^^ + dA^ - 2a2D - 2mc2F* = 2{m - D){1 - D)B^ , (197) 

- dDA^ - mh2DF^ + aal)^ + m^C2F'' = D{1 - D)B'^ . (198) 

On solving these equations we obtain the four unknowns A,F,B,D and further also obtain two constraints 
between the seven parameters. We find 

2 2B'^[b2{2mD -m + 2D- ZD^) - d{m + 2D- 2mD - D'^)] 



A' 



(6162 - d2) 



2 _ 2B^[bii2D + m - 2Dm - D^) - d{2mD - m + 2D - 3D^)] 

{m - D){b,b2 - d^) ' 

y 3D(1-D)(m-D) , 9 r 1 9 . X 

5 = rr-, ^ ^ , oi + dF^ = 130-1- m]B^ , 200) 

4ci[62(2mL>-m + 2Z)-3L>2)-d(m + 2D-2mD-L»2)]2 ' ^ J > v ; 

while the two constraints are 

ci^^ = {m- Dfc2F^ , (201) 

, 5 B'^Um - Df + 3D(1 - D)] , , 

02 - 62F' + / = . 202 



27 



Prom positivity considerations, one can show from here that 



3 

Further, it fohows that ai < if < (1 + m)/3. 

At m=l, the solution H192|) goes over to the hyperboUc, bright-dark sohton solution H181|) and the 
constraints (|199j) to (|202|) take the simpler form as given by Eqs. (|182j) to (|185j) . 

Special case of 6162 = (P 

One can show that the solution 1)1921) exists even in the case 6162 = d?. It turns out that such a solution 
exists only if 

V^[2(l + m)D -m- SD^] = y^[m - D"^ + 2(1 - m)D] , (204) 

and further while relations H2U0|) to ()202|) are still valid, the relation 1)199(1 is no longer valid. Instead of 
separate expressions for A^, F"^, one now only has a constraint 



biA^ + VM2(1 - D)F^ = 2B^[2{1 + m)D - 3D^ - m] . (205) 

In the limit m ^ 1, the hyperbolic bright-dark solution ()181|) is still valid provided the constraints 
(|TH5|) and (fTHU]) as well as Eqs. (fTHH|) to ((TH5|) are satisfied. 

Energy: Corresponding to the periodic solution (|192() . the energy E and the constant C are given by 



1 



C 



2 3 2 3 



(206) 



2D^ 

where /i,/3 are as given by Eqs. ()103() . I|1U4|) . (|152|) and (|153|) . On using the expansion formulas for 
E{k), K{k) and n(D, k) around = 1 as derived above, the energy of the periodic solution can be rewritten 
as the energy of the corresponding hyperbolic (bright-dark) soliton solution [Eq. (|181|) ] plus the interaction 
energy. We find 

E = Ek^nk + E,nt = BA^[l['^ + A;'2/«] + SF^i/f + A;'^/^] , (207) 

where are given by Eqs. (fTT^ and (|Tni) while if^ = and are given by Eqs. H136|) and H155|) . 
respectively. 

Note that this solution exists only when &162 > (P. The interaction energy vanishes at /c = 1. 
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2.13 Solution XVII 

So far, we have presented ten solutions at T = Tc and six solutions for T < Tq- We now present one 
solution which exists when Tc < T < Tp where Tp denotes the point of inflection where cp^ potential (at 
d = 0) has an absolute minimum at (/» = and two points of inflection. Note that the point of inflection 
occurs when 6^ = 4aiCi. 
It is easily shown that 

A F 

, ^ , (208) 



•yr^^"Dsn^^^Bx"+~r()~m^ ' ■yr^^TDsn^^SaT+'xoyrTi^ 
is a solution to the field Eqs. @ provided the following six equations are satisfied 

ai - biA^ + ciA^ + dF^ = DB^ , (209) 

- dF^ + biA^ - 2ai = 2{D -I- m)B^ , (210) 

aiD = {3m-D- Dm)B'^ , (211) 
02 - 62^2 + C2F^ + dA"^ = DB^ , (212) 

- dA'^ + b2F^ - 2a2 = 2{D -I- m)B^ , (213) 

02-0 = {3m-D- Dm)B^ . (214) 

Prom here one can determine the four unknowns A, F, B, D and further one obtains two constraints between 
the seven parameters 01^2, ^i,2i ci^2j d. We get 

^2 aiD 2 2{d + b2)B^[D^ -2{l + m)D + 3m] 

= T -pr-: TT , A 



[3m-D{l + m)]' D{bib2-d'^) 



,2 2{d + bi)B'^[D'^ -2(1 + m)D + 3m] 



D{bib2 - d2) 
[1)2 _ 2(1 +m)Z) + 3m]2 8(6162 - d^)^ 



(215) 



(1 - D){m - D)[3m - D{1 + m)] 4.aici{d + 62)2 ' 
while the two constraints are 

01=02, C2(d + 6i)2 = ci(d + 62)^ (216) 
Note that 01,02 are always positive. 
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In the limit d = 0, one obtains the uncoupled pulse lattice solution PH [T^ [T^ 

A 



satisfying 



a/I — Dsii^{Bx + xo,m) 



^2 252[D2-2D(l+m) + 3m] {I - D)ai 

A = — , B — 



(217) 



Dbi ' [2m - 1 - D{2 - m)\ 

[1)2 _ 2D{1 + m) + 3m]2 36f 



(l-D)(m-£))[3m-(l + m)Z)] 4aiCi ' 
At m = 1 the solution reduces to the nontopological (dark-dark) solution 



(218) 



A F 

, V= , (219) 



provided 



1- Dtanh\Bx + xo) Jl - D tanh^Bx + xq) 



(3 - 2D) ' D{bib2 - (P) 

n)B^il-D){3-D) (3-D)' _ 
D{bib2-cP) ' (3-21?) 4aiCi(d + 62)2 

while the two constraints are the same as given by Eq. H216() . The constraint < Z) < 1 implies the 



^2 2{d + b^)B^{l-D){'i-D) ^ (3-D)2 2>{b,b2-(f? ^ ^220) 



relation 

4aici(d + 62)^ < (&162 - d2)2 < (16/3)aici(d + 62)^ , (221) 
which for d = reduces to the constraint 

4aici < 6? < (16/3)aici , (222) 

i.e. Tc < T < Tp. This corresponds to the situation when there is an absolute minimum at (j) = and there 
are two degenerate local minima |111 1121 [T^ . 
Special case of 6162 = d? 

One can show that the solution (|2U8|1 exists even in the case 6162 = d? provided m < 1. It turns out 
that such a solution exists only if 

biA^ = b2F^ , C2bi = C162 , (223) 
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and further 



D = {1 + m) — \/l — m + , ciA^ 



3^2(1 - D){m- D) 



B^[3m - D{1 + m)] 



(224) 



Energy: Corresponding to the periodic solution (|2U8() . the energy E and the constant C are given by 



(using appropriate integrals in ^7j) 



where 



with 



m 



C = ^{A^ + F^)B\ 



[GiK{k) + G2E{k) + G^Yl{D, k)] 
4D^{1- D){k'^ - D) 



(225) 



(226) 



Gi = D{k'^ - D){D'^ -2D + ADk'^ - 3k'^) , G2 = D'^[D'^ - 2D{1 + A;^) + 3k'^] , 

G3 = 4D{1 - D){k^ - D)[3k^ - (1 + k^)D] - D[D^ - 2(1 + k^)D + 3k^f . (227) 

On using the expansion formulas for E{k), K{k) and 11(1), k) around k = 1 derived above, the energy 
of the periodic solution can be rewritten as the energy of the corresponding hyperbolic (dark-dark) soliton 
solution [Eq. (|219|) ] plus the interaction energy. We find 



E = Ek^nk + Eint = (A^ + F^)B[ir + k'^l 



4 J ' 



(228) 



where I^^^ and I^^ are given by 



r(0) 



(3 -D) {D + 3) 



4(1 - D) 8L>/2 



In 



1 + 

1-VD 



^(1) _ -(Zj2_8D + 3) , 3{D + 4) 

161)3/2 



16(1 -Z))2 



+ 



In 



1 + VD 



1 



(229) 
(230) 



Note that this solution exists only when > d^. The interaction energy vanishes at = 1. 
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2.14 Solution XVIII 



Finally, we present three novel mixed phase solutions, i.e. in which one of the field exists for T > while 
the other one exists for T < Tc. 
It is easily shown that 

A FsniBx + xo,m) ^^Sl) 



v^l — Dsn^{Bx + xq, m) ' ^/l — Dsn^{Bx + xq, m) 

is a solution to the field Eqs. © provided the following six equations are satisfied 

ai - biA^ + ciA^ = DB^ , (232) 

dF^ + biDA'^ - 2aiD = 2D{D - 1 - m)B'^ , (233) 

aiD - dF"^ = {3m- D- Dm)B'^ , (234) 

02 + dA^ = (3L» - 1 - m)B'^ , (235) 

- dDA^ - b2F^ - 2a2-D = 2{m -D- mD)B^ , (236) 

a2D^ + b2DF^ + caF^ = mDB^ . (237) 

On solving these equations, we obtain the four unknowns A, F, B, D as well as two constraints between the 
seven parameters. We find 

A^ 



2 2B'^[b2{D'^ -2D- 2mD + 3m) + d{2D + 2mD -m- SD"^)] 



{bib2 - d^) 

'?j^'^\d(n'^ -9.n -2m n 

F 



2 2B'^[d{D'^ -2D- 2mD + 3m) + bi{2D + 2mi:> - m - SD"^)] 



{bib2 - d^) 

^2^ 3{l - D)im - D)ibib2 - d^ 

ADci[b2{D^ -2D- 2mD + 3m) + d{2D + 2mD - m - 3D^)]^ ' 

aiD - dF'^ = [3m - (1 + m)D]B'^ , (238) 
while the two constraints are 

D^cU^ = C2F^ , 02 + dA^ = (3D - 1 - m)B'^ . (239) 

At m = 1 the solution reduces to the dark-bright solution 

A Ft.nHBx + xo) ^^^^^ 



1- Dtaiih''{Bx + xo) ^1 - tanh^(5x + xq) 
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provided 

2 _ 2^^(1 - D)[b2{3 -D) + d{3D - 1)] 
(6162 

2 ^ 2^2(1 - D)[d{3 -D) + h{3D - 1)] 
(6162 - d2) 

- .nc.lj%fdln-iw ' - ^ - ^^^^^ ' ^'''^ 

while the two constraints are 

D^ciA'^ = C2F^ , 02 + dA^ = {3D - 2)B^ . (242) 

Energy: Corresponding to the periodic solution ()231() . the energy E and the constant C are given by 
(using appropriate integrals in \17\ ) 

E = {A^h + F^h)B , 

C = !^A'B' - ^[B\l + m-D)+a2], (243) 

where h and h are given by Eqs. (fTUH|) . (|TU11) . and (P77|) . 

On using the expansion formulas for E{k), K{k) and Il{D, k) around /c = 1 as derived above, the energy 
of the periodic solution can be rewritten as the energy of the corresponding hyperbolic (dark-bright) soliton 
solution [Eq. H240|l ] plus the interaction energy. We find 

E = Ek^nk + E,^t = BA^[lf^ + + BF^[lf^ + A;'^/^] , (244) 

where are given by Eqs. (ITT^ . (tM . and (0^171) . 

Note that this solution exists only when 6162 > d'^ . The interaction energy vanishes at k = 1. 

2.15 Solution XIX 

Another novel mixed phase solution is given by 

^ ^ ^ ^ ^ Fcn{Bx + XQ,m) 

y^l- Dsn^{Bx + xo,m) v/1 - Dsn^{Bx + xo,m) 

This is a solution to the field Eqs. Q provided the following six equations are satisfied 

ai - biA^ + ci^^ + dF^ = DB'^ , (246) 
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- d{l + D)F^ + biDA^ - 2aiD = 2D{D - 1 - m)B^ , (247) 

aiD + dF"^ = (3m-D- Dm)B'^ , (248) 

a2 + dA'^-b2F^ + C2F'^ = -{l-D)B^, (249) 

- dDA^ + 62(1 + D)F'^ - 2c2F^ - 2a2-D = 2(1 - D){m - D)B^ , (250) 

a2-D^ - b2DF'^ + C2F^ = mD{l - D)B'^ . (251) 

On solving these equations, we obtain the four unknowns A, F, B, D as well as two constraints between 
the seven parameters. We find 

2 _ 252[62(L>2 -2D- 2mD + 3m) - d{m + 2mD -2D- D"^)] 

F^ = 



B^ = 





D{bib2 - d^) 




2B'^[bi{m + 2mD - 


2D - D^) - d{D^ - 


2D - 2mD + 3m)] 




(l-L>)(M2-d2) 




3D{1 


-D){m-D){bib2- 


-d^f 



4ci[62(£>2 _ 2D - 2mD + 3m) - d{m + 2mD -2D- L»2)]2 ' 
aiD + dF"^ = [3m - (1 + m)D]B'^ , (252) 

while the two constraints are 

9 .4 . x9 d 9 B^\D{mD + X)-2m(l + D)] , , 

D^cxA^ = (\-Dfc2F\ a2D-b2F^ = ^^ _ ^ ^. (253) 

At m = 1 the solution reduces to the dark-dark soliton solution 

, A , Fsech(^a; + xp) 

^= I ^= I (254) 

\Jl-D tanh^(Sx + xq) ^Jl - D tanh^(Sx + xq) 



provided 



2 _ 2B^{1 - D)[b2{3 -D)- d{3D - 1)] 
£'(6162 - d^} 
^^^2B^Ml + D)-d{3-D)] 



(6162 - (P) 

while the two constraints are 

D'^ciA'^ = (1 - Dfc2F'^ , a2D - 62-^^ + (2 + D)B'^ = . (256) 
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Special case of 6162 = 

One can show that the solution (|245|1 exists even in the case 6162 = d?'- It turns out that such a solution 
exists only if 

V^[L>2 _ 2(1 + rn)D + 3m] = V^[m(l + 2D) - D{2 + D)] . (257) 

At m = 1 this implies 

^6^(3 -D) = + D) . (258) 

Energy: Corresponding to the periodic solution (|245() . the energy E and the constant C are given by 
(using appropriate integrals in [T7] ) 



"^2(1 - mD) 



C = —A^B^ - — 



{l-D) 



02 



(259) 



2 4.D 

where h and I4 are given by Eqs. (fnH|) . (ITO . (1^ and (1^ . 

On using the expansion formulas for E{k), K{k) and Il{D, k) around A; = 1 as derived above, the energy 
of the periodic solution can be rewritten as the energy of the corresponding hyperbolic (dark-dark) soliton 
solution [Eq. (|254|l ] plus the interaction energy. We find 

E = Ek,nk + Eint = BA^[lf^ + k'^lf^] + BF^lf^ + fc'^/f )] , (260) 

where 4°/^ are given by Eqs. (tTHHl) . (tTTTl) . (tCT and (OHm . 

Note that this solution exists only when 6162 > d'^- The interaction energy vanishes at k = 1. 

2.16 Solution XX 

Finally, yet another novel mixed phase solution is given by 

A F(ln{Bx + XQ,m) 

, ■0 = ^ , • (261) 



^1 — Dsn^{Bx + xq, m) ' ^/l — Dsn-^ (i? + xq , m, 

This is a solution to the field Eqs. @ provided the following six equations are satisfied 

ai - 61^2 + ciA^ + dF^ = DB^ , (262) 

- d{m + D)F'^ + hiDA^ - 2aiD = 2D{D - 1 - m)B^ , (263) 
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aiD + dmF'^ = (3m -D - Dm)B'^ , (264) 

a2 + dA^ -b2F'^ + C2F^ = -{m- D)B^ , (265) 

- dDA^ + b2{m + D)F'^ - 2mc2F* - 2a2D = 2(1 - D){m- D)B'^ , (266) 

a2D^ - b2mDF^ + C2m?F^ = D{m - D)B^ . (267) 

On solving these equations, we obtain the four unknowns A, F, B, D as well as two constraints between 
the seven parameters. We find 

,2 _ 2B^[b2{D'^ -2D- 2mD + 3m) - d{m + 2D- 2mD - D^)] 

D{bib2 - d^) ' 
2 _ 2B'^[bi{m + 2D- 2mD - D"^) - d{D'^ -2D- 2mD + 3m)] 

{m - D){bib2 - d^) ' 

2 3D{1 - D){m - D){bib2 - d^f 

(Z)2 -2D- 2mD + 3m) - d{m + 2D- 2mD - D'^)f ' 

aiD + mdF'^ = [3m - (1 + m)D]B'^ , (268) 



A 
F 
B 



while the two constraints are 

^9 .4 . ^4 , ^9 B'^ \ D"^ + mD( 3m - 2) - 2771^] 

D^cxA^ = (m - Dfc2F^ , a2D - mbiF^ = — ^ ^ ^ . (269) 

B 

At m = 1 the solution reduces to the dark-dark soliton solution (|254|) satisfying the relations H255j) and 
(EHHl)- 

Special case of 6162 = d"^ 

One can show that the solution 1)2611) exists even in the case 6162 = d'^. It turns out that such a solution 
exists only if 

VhiD^ - 2(1 + m)D + 3m] = ^/h[m + 2D - 2mD - D^] . (270) 

At m = 1 this implies the constraint (|258j) . 

Energy: Corresponding to the periodic solution ()261j) . the energy E and the constant C are given by 
(using appropriate integrals in |17j l 



E = {A^h + F^l3)B 

C = —A^B^ 

2 4D 



mB'^(2m — vn? — D) 
{^^^ "2 
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(271) 



where I3 and I4 are given by Eqs. H152|) and 1)2261) . respectively. 

On using the expansion formulas for E{k), K{k) and 11(1), k) around k = 1 derived above, the energy 
of the periodic solution can be rewritten as the energy of the corresponding hyperbolic (dark-dark) soliton 
solution [Eq. (|254p ] plus the interaction energy. We find 

E = Ek,nk + Eint = BA^[li'^ + A;'2/f )] + BF'[4^^ + k'^I^^^] , (272) 

where /f^ = 4°\4^\/f '^^ are given by Eqs. ([THH|) . (IT^ . XTM and 

Note that this solution exists only when 6162 > d?- The interaction energy vanishes at A; = 1, as it 
should. 

Interaction energy: Summarizing, we thus have obtained twenty solutions of the coupled (p^ model. 
We have calculated the corresponding soliton and interaction energy in each case. If we carefully look 
at the interaction energy expression in each case, we find that except for the solution XI, for all other 
solutions, the soliton interaction energy (in the dilute gas or asymptotic limit J3]) has the leading term 
of the form k''^ ln{A/k'). Out of these nineteen solutions, solutions I, II, III, XII and XV have repulsive 
interaction energy while the solutions IV, V, VI , XIII, XVI and XIX have attractive interaction energy. 
However, for the solutions VII, VIII, IX, X, XIV and XX, the interaction energy can be attractive or 
repulsive, depending on the values of the various parameters of the model. Finally, for solution XI, the 
interaction energy (in the dilute gas or asymptotic limit) has the leading term of the form /c'^ and is 
repulsive. Since the lattice periodicity r = 2nK{k)/B with n = 1 for dn based solutions, and n = 2 for sn 
and cn based solutions, using K{k) ~ ln(4/A;') we find that k'^ = 4n^ exp(— ri?/n). Thus the interaction is 
purely exponential for solution XI whereas it is of the form fe'^ ln(4:/k') = (2rnB) eyip{—rB/n) for all other 
solutions. 

3 Solutions of a Coupled 0^ Discrete Model 

We shall now consider a discrete variant of the above continuum coupled (p^ model and obtain its static 
solutions. We start from the coupled static field Eqs. @. The discrete analogue of these field equations 
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has the form 

{l/h?)[4)n+l + 4>n-l - '^4>n\ = ai(pn " ^1'/'^ + Cl</>^ + dlp^Cpn, 
{l/h^)[tpn+l + V'n-l - 2-071] = a2'lpn " b2V'n + ^l^^n + d(t)n'>Pn, 

(273) 

where /i is the lattice spacing. We are unable to find any solution of these coupled equations. However, as 
in the Ablowitz-Ladik discretization of the nonlinear Schrodinger equation if we replace 4),ip in the 
terms by the corresponding average, then we can find the exact solutions. Thus, we consider the 
coupled equations 

[(/>„+! + - 24>n] = ai0„ - hi4>l^ + -J<Pn[(Pn+l + 4>n-l] + d'lpl<j3n , 

{l/h'^)[lpn+l + tpn-1 - 2V'n] = a20n " &20i^ + yV'n[V'n+l + V'n-l] + ^^V^n • 

(274) 

We shall show that while in the continuum case there are twenty periodic solutions, we are able to obtain 
only six periodic solutions in the discrete case. We shall also see that none of the discrete solutions has 
a smooth continuum limit, i.e. in the limit h ^ 0, none of them goes over to the continuum solutions 
obtained in the last section. 

3.1 Solution I 

It is easy to show that 

(pn = Ay/rnsn{hB[n + c],m) , ipn = F^/rnsn{hB[n + c],m) , (275) 

is a solution to the above equations provided 

ai = 02 = (2//i2) [cn(/i5, m)dn{hB, m) - 1] < , (d + bi)A'^ = {d + 62)-?^^ , 
2 2((i + b2)sn^{hB, m)cn{hB, m)dn{hB, m) 



c,A^ = c,F^ , era, = [--ihB,m)MhB m) - l]{d- - h,h,f 

cn^{hB,m)dn^{hB,m){b2+d)'^ ^ ' 
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It is worth noting the completely different form of the solutions in the continuum and the discrete cases. 
We shall see that this is true for all the six solutions that we obtain in the discrete case. Also note that 
while in the continuum case 6162 > d^, in the discrete case cP > 6162. 
m = 1 limit 

In the 171 = 1 limit, this solution goes over to the bright-bright topological solution 



provided 



Atanh(/i5[n + c]), ^p = F tanh{hB[n + c]) , (277) 



a, = a, = -i2/h^) t.nl.HhB) < , = ^+ tanh^(/.5)sech^(/.i?) 



t&n\\^{hB){(P - bib 



^2 



ciai = ^ , (278) 

(52 + (i)2sech^(/iS,m) ^ ' 

while the other two relations are as given by Eq. H276() . We note that this solution exists only if dP > 6162- 
In the limit d = (i.e. uncoupled case), the solution is given by 

(t)n = Ay/msn{hB [n + c] , m) , (279) 

provided 

ai = {2/h^)[cn{hB,m)dn{hB,m) -1] <0, = ^ ^ ' v v , j 



h?bi 



ciOi [cn(/ii?, m)dn(/ii?, m) — 1] 
b\ crP {hB,m)diP{hB,m) 

In the limit m = 1 the uncoupled solution is given by 



(280) 



= ^tanh(/i5[n + c]), (281) 

provided 

_ 2tanh2(/iB) 2 _ 2tanlP{hB)sech^{hB) ciUi _ tanlP{hB) 

JP ' ^ " hH[ ' 'W~~^ech\hB) • ^^^^^ 

Thus, such a solution exists provided ai, 61 < 0. 
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3.2 Solution II 



We shall now obtain three periodic solutions, all of which in the limit m — > 1 go over to the same dark-dark 
soliton solution. The first one of these, i.e. 



A^/rncn(hB [n + c],m) , ip = F^/mcn{hB [n + c] , m) , 



(283) 



is a solution to the above equations provided 

cn{hB, m) 



oi = a2 = (2//i 



dn^(/i5,m) 



2 2{d + b2)sn^{hB,ni)cn{hB,m) 



A 



h?{hih2 -cP)dn*{hB,m) 



ciai 



cn{hB, m) 
dn^ {hB,m) 



1 



dn4(/iS,m)(6i62 -d- 



2\2 



(284) 



cn'^{hB,m){b2 + dy ' 

Again notice the completely different form of the solutions in the continuum and the discrete cases. One 
can show that ai, 02 > (<) provided m > (<) 1/2. Note that this solution exists only if 6162 > d'^. 

In the m = 1 limit, three of the solutions (i.e. solutions II, as well as solutions III and IV to be discussed 
below) go over to the nontopological soliton solution 



Asech{hB[n + c]) , ipn = Fsec\i{hB[n + c]) , 



(285) 



provided 



ai = 02 = {2/h^)[cosh{hB) - 1] > , = 

id + bi)A^ = id + b2)F\ ciA^ = C2F\ 

, , „ ,sech'^(hB)(bib2-d'^)^ 
ciai = [cosh(/iB) - 1]— ^ ^ 



2 2{d + b2)smh^{hB)cosh{hB) 



h?{d?-bib2) 



{h2 + dY 

In the limit d = (i.e. the uncoupled case), the solution is given by 



(286) 



4>n = A\/mcn{hB[n + c], m) , 



(287) 



provided 



ai = (2//i2 



cn{hB, m) 
dn^ {hB,m) 



A 



2 2sn'^{hB,m)cn{hB,m) 



h'^bidn'{hB,m) 



cn(hB, m) 
dn^ {hB,m) 



- 1 



dn^ {hB,m) 
cn^ihB^m) 



(288) 
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In the limit m = 1, the uncoupled solution takes the form 



Asedi{hB[n + c]) , 



(289) 



provided 



ai = (2//i2) [cosh{hB) - 1] , 



2tanh2(/iS) 
h%isech^{hB) ' 



([1 -sech(/i5)]sech(/iB)) . 



Thus, the uncoupled solution exists only if ai, 6i > 0. 



3.3 Solution III 



It is easy to show that another periodic solution is given by 



(290) 



Adii{hB[n + c], m) , ip = Fdn{hF[n + c], m) , 



(291) 



provided 



dn{hB, m) 



{d + bi)A^ = {d + b2)F^ , ciA^ = C2F^ 



>0, A 



2 2{d + b2)sn^{hB,m)dn{hB,m) 



h^{bib2-d^)cn^ihB,m) 



ClOl 



dn{hB, m) 
cn^ {hB, m) 



1 



'2\2 



cii^{hB ,m){bib2 — d 
drP'{hB,m)(b2 + dY 



(292) 



Note that as in the continuum case, here too ai,02 > 0. Further, this solution exists only if 6162 > ■ 
Finally, in the limit m = 1, this solution goes over to the solution H285|) satisfying the constraints (|286() . 
In the limit d = (i.e. the uncoupled case), the solution is given by 



Adn{hB [n + c] , m) 



(293) 



provided 



dn{hB, m) 
cn^{hB, m) 



A 



2 2sn'^{hB,m)dn{hB,m) 



h%icn'^{hB,m) 



dn{hB, m) 
cn^{hB, m) 



1 



c\i^{hB, m) 
dn^{hB,m) 



(294) 



In the limit m = 1, the uncoupled solution goes over to the solution H289() satisfying relations H29U() . 
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3.4 Solution IV 



It is easy to show that yet another periodic solution is given by 



A^/rncn(hB[n + c],m) , ip = Fdn{hB[n + c],m) , 



(295) 



provided 



cn{hB, m) 
dn^ {hB,m) 



,a2 = (2//i2) 



dn(hB, m) 



cn^ (hB, m) 



>0, 



A 



2 _ 2[ddn^{hB,m)+b2cn^{hB,m)]sn^{hB,m) 
/i2(6i62 - d2)cn4(/i5, m)dn^(/iB, m) 
2[6idn^(/i-B, m) + dcn^{hB, m)]sn'^{hB, m) 



h?{hih2 - d^)cn^{hB,m)dn'^{hB,m) 

C2F' 



4 2sn'^{hB,m) ^4 2sn'^{hB,m) 
c\A " 



(296) 



h?dn^{hB,m)' h?cn^{hB,m)' 
Note that 02 > while ai > (<) depending on if m > (<) 1/2. 

Again note that like the previous two solutions, this one also exists only if 6162 > <P . As mentioned 
above, in the limit m = 1, this solution also goes over to the solution H285|) satisfying the constraints H286() . 
Further, in the decoupled case (i.e. d = 0), (j)n and 'i/'n satisfy solutions (|287|1 and H293|l . respectively, with 
the appropriate constraints. 



3.5 Solution V 

Finally, there are two periodic solutions, both of which at m = 1 go over to the (same) bright-dark soliton 
solution. The first one is given by 



Ay/rrisTi.{hB[n + c], m) , -0 = F\/rncn{hB[n + c],m) , 



provided d < and further 



(297) 



ai - 'm\d\F'^ = {2 / h^)[cn{hB , m)dn{hB, m) - 1] < , 



a2-m\d\A^ = {2/h^ 



cii{hB, m) 
dn^{hB,m) 



A 



2 2[\d\ - b2dn^{hB, m)]sn'^{hB, m)cn{hB, m) 



h'^{hib2-d^)dn'{hB,m) 
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^2 _ 2[6i - \d\dn''{hp,m)]sn^{hp,m)cn{hp,m) 
~ /i2(6i62 - d2)dn^(/iS,m) ' 

4 2sn^(/i.B,m) 4 2sn^(feB,m) 

n2 h^dn {hB,m) 

Note that this solution exists only if 6162 > c^^- In the decoupled limit (i.e. d = 0), the (pn and ^/^^ fields 

go over to the solutions (|279jl and l|287|l . respectively, satisfying appropriate constraints. 

3.6 Solution VI 

Finally, the second solution is given by 

(j) = A^/msn{hB[n + c\,m) , ijj = Fdii{hB[n + c],m) , (299) 
provided d < and further 

ai - |d|F2 = {2/h'^)[cn{hB,m)dn{hB,m) - 1] < , 
02 - \d\A^ = {2/h^) 



dn{hB,m) ^ 



> 0. 



A 
F 



cn'^{hB, m) 

2 2[|d| - b2cn^{hB,m)]sn^{hB,m)dn{hB,m) 

/i2(6i62 -d2)cn4(/i5,m) 
2 _ 2[6i - \d\cn^{hB,m)]sn'^{hB,m)dn{hB,m) 



h?{bih2 - d^)cn^{hB, 



m] 



.A 2sn'^(hB,m) . 2sn^(hB,m) 

= ^2F = /.2 4\.R V 300 

n^cn*(n,i>, mj 

Note that 02 > 0. 

In the limit m = 1, both of these solutions (i.e. solutions V and VI) go over to the bright-dark soliton 
solution 

^ = Atanh(/iB[n + c]) , tp = Fsech{hB[n + c\) , (301) 

provided d < and further 

ai - \d\F^ = -(2//i2)tanh2(/iS),a2 - \d\A'^ = {2 / h^)[cosh{hB) - 1], 
2 2[|d| - 62sech5(/i5)] tanh2(/iS) 



/i2(6i62 -(i2)sech3(/iS,m) ' 



,2 _ 2[bi - \d\sech^ihB)]tanh^{hB) 



h?{hib2 - d'^)sec\i\hB) 
/i2 ' "^2^ /l2 



^^^4^2tanh^^ ^^^,^2sinHhB)^ ^^^^^ 
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Note that this solution exists only if 6162 > d^- In the decoupled limit (i.e. d = 0), the and the V'n 
fields go over to the solutions (|279|) and l|293|) . respectively, satisfying appropriate constraints. 

As in the continuum case, can one say if these discrete solutions are valid below, above or at Tc? There 
is a note of caution here. Since we are considering a discrete model here, and since the corresponding 
Hamiltonian from which the coupled discrete field Eqs. (|274|) can be derived is not known, strictly speaking 
we cannot directly associate oi, 61 and other parameters of the discrete problem and draw conclusions about 
the nature of the transition by using the corresponding continuum model potential 

However, as a first guess, it might be reasonable to assume that for the same values of the ratio of the 
parameters (i.e. aici/bl) as in the continuum case, one can classify if the solution is below, above or at 
Tc- For example, since solution I exists if ai,6i < (ci is always assumed to be positive), this implies 
that this solution exists with T < T/^. It is worth reminding that in case 61 < 0, then the uncoupled (j)^ 
model characterized by the potential corresponds to a second order (and not a first order) transition. 
Similarly, since for solutions II, III and IV, 61, ai > and > 4aiCi, this implies that T < T^. Finally, 
solutions V and VI correspond to mixed-type solutions with field corresponds to T < T/^ while ipn field 
corresponds to T < . 

Since we do not know the Hamiltonian associated with Eqs. 1)274(1 . we are unable to find the soliton and 
interaction energy explicitly for these discrete solutions and hence cannot comment whether the Peierls- 
Nabarro barrier is zero or non-zero. 

4 Conclusion 

We have provided a set of twenty distinct exact, periodic domain wall solutions for a coupled (j)^ model with 
biquadratic coupling. These in turn lead to nine distinct (hyperbolic) soliton solutions both at, above and 
below Tc- The corresponding discrete case was also considered and we found six different periodic solutions 
which in turn lead to three distinct (hyperbolic) soliton solutions. The soliton and interaction energy were 
obtained for all twenty continuum solutions. For the six solutions of the discrete model, the calculation of 
the Peierls-Nabarro barrier j2Ul 1211 122j and soliton scattering |231 are important topics of further study. 
Similarly, scattering of solitons for the twenty solutions in the coupled (j)^ continuum model is an interesting 
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issue with these static solutions boosted with a certain velocity. We also discovered previously unknown 
periodic solutions of the uncoupled (j)^ model [TTl [T^ [T^ . 

It would be useful to explore whether the different solutions are completely disjoint or if there are 
any possible bifurcations linking them via, for instance, analytical continuation. We have not carried 
out an explicit stability analysis of various periodic solutions P^. However, the energy calculations and 
interaction energy between solitons (for m ~ 1) could provide useful insight in this direction. Similarly, 
it would be worth exploring the problem of a coupled (j)^ model in the presence of an external field. Our 
results are relevant for domain walls in structural phase transitions in ferroelectrics and elastic materials 
[HI El ^1 ^1 ^1 1^ and possibly in field theoretic contexts [HI d El- These ideas and solutions can be 
generalized to other coupled models for first order transitions such as coupled asymmetric double well 
potentials (cp'^-cp^-cp'^) and will be discussed elsewhere. 
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